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Abstract 


Work  in  the  period  under  review  centered  principally  on 
the  topics:  (1)  numerical  methods  for  systems  of  linear 
equations,  (2)  the  numerical  range,  including  computer 
routines  for  its  evaluation,  (3)  inequalities  involving 
matrices,  f 4 )  properties  of  eigenvalues,  sinqular  values,  and 
invariant  factors,  (5)  properties  of  norms  and  condition 
numbers  on  matrices. 

Under  (1) ,  p-adic  arithmetic  (from  pure  number  theory) 
has  now  been  found  to  be  well  suited  for  computing  exact 
solutions  of  linear  systems  that  may  be  ill  conditioned,  and 
work  is  under  way  to  develop  computer  routines  to  implement 
this  new  procedure.  For  (2) ,  a  blend  of  analytical  and 
computer  techniques,  supplemented  by  geometrical, 
considerations,  continues  to  yield  new  insights  concerning  the 
various  numerical  ranges  and  numerical  radii  on  matrices. 
This  research  is  motivated  by  applications  to  partial 
differential  equations.  (3)  Insightful  understanding  of 
numerical  aspects  of  linear  algebra  requires  insightful 
understanding  of  the  numerous  inequalities  that  pervade  the 
subject.  Continued  progress  in  this  direction  of  research  is 
reported,  there  being  inequalities  involving  permanents, 
nonnegative  matrices,  submatrices  of  Hermitian  matrices, 
diagonal  elements  of  normal  matrices,  matrix  valued 
inequalities,  and  others.  (4)  Eigenvalues,  singular  values. 


invariant  factors,  important  in  every  application  of  .linear 
algebra,  continued  under  investigation,  with  new  results 
reported  for  matrices  of  algebraic  integers,  and  for  normal 
matrices.  (5)  Norms  and  condition  numbers,  central  tools  for 
the  numerical  analyst,  remain  under  active  investigation,  with 
new  results  connecting  von  Neumann  norms  to  condition  numbers, 
and  also  new  results  for  the  very  intransigent  multi- 
plicativity  factors  that  make  norms  submultiplicative. 
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ABSTRACT 


MARVIN  MARCUS:  MATRIX  NORMS,  CONDITION  NUMBERS ,  AND  THE 
NUMERICAL  SOLUTION  OP  LARGE,  SPARSE,  LINEAR  SYSTEMS 


The  author's  research  during  the  period  October  1,  1981 
through  September  30,  1982  was  primarily  concerned  with  the 
following  topics:  fl)  inequalities  for  unitarily  invariant 
norms  and  condition  numbers;  (2)  convexity  properties  of 
certain  eigenvalue  inclusion  regions  arising  from  quadratic 
and  bilinear  forms;  (3)  inequalities  for  principal  submatrices 
of  positive  definite  hermitian  matrices;  f 4 )  generalized 
matrix  functions  and  the  index  of  a  svmmetry  class  of  tensors. 
Activities  included  the  publication  of  three  papers  in  these 
areas  and  the  presentation  of  some  of  these  results  at  the 
April  1982  conference  on  Applied  Linear  Algebra  held  by  the 
Society  for  Industrial  and  Applied  Mathematics  (SIAM) . 


The  present  annual  report  summarizes  this  research  of  M. 
Marcus  sponsored  by  the  Air  Force  Office  of  Scientific 
Research. 


cm* caw 
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1.  Inequalities  for  Uni tarily  Invariant  Matrix  Norms  and 
Condi tion  Numbers  . 

It  is  pointed  out  in  the  LINPACK  User's  Guide  [31  that 
obtaining  estimates  of  norms  and  condition  numbers  is  decisive 
in  dealing  with  two  general  problems  that  often  arise  in  the 
numerical  solution  of  large  linear  systems. 

If  A  is  an  n-square  complex  matrix,  how  close  is  A 
to  being  singular? 

How  much  does  the  solution  change  when  A  is  reolaced 
by  A+E,  where  E  is  a  perterbation  matrix 
representing  the  total  error  arising  from  initial 
data  errors,  round-  off  errors,  etc.? 

A  well  known  theorem  r.121  states  that  the  relative  error  in 
A+E  may  be  magnified  bv  as  much  as  the  condition  number  of  A, 
c  (A) .  Recall  that  if  JjAj]  is  a  matrix  norm  then  the  condition 
number  is  defined  by 

C(A)  -  ||  A  |[  ||  A"1 1) 

In  order  to  obtain  computable  lower  bounds  on  c(A),  the 
investigator  recently  proved  that  for  the  Hadamard  product 


A-B  -  [a44b..  J 
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the  following  inequality  holds 

II  A-B  II  £  ||a||  || B  j| 

where  the  indicated  norm  is  the  usual  matrix  norm  subordinate 
to  the  Euclidean  vector  norm.  We  thus  have  that 

II  A-  A_1|l  £  c(A).  (1) 

We  have  extended  this  result  to  a  large  class  of  von  Neumann 
norms  [10].  For  example,  if  A  is  an  n-square  matrix  with 
singular  values 

.  a.  st  •••  it  a  st  o 
l  n 

and  if 

p  k  1/p 

II  All  =  Jj  A||  S(ZaMP 
k  m=l  m 

then  fl)  holds.  Research  is  still  underway  to  determine 
whether  or  not  this  result  holds  for  all  unitariltv  invariant 
norms  and  in  what  cases  equality  can  hold  for  (l). 


Arising  from  Quadratic  and  Bilinear  forms  . 

Let  A  be  an  n-square  complex  matrix.  The  numerical  range 
of  A,  W (A) ,  is  defined  as  the  image  of  the  unit  sphere  under 
the  quadratic  mapping 

x  — >  x*Ax. 

A  classical  result  due  originally  to  F.  Hausdorff  and  0. 
Toeplitz  F4,161  states  that  the  numerical  range  of  A  is  always 
a  convex  region  in  the  plane  which  contains  the  spectrum  of  A. 
Thus  W (A)  can  be  used  to  yield  estimates  of  the  eigenvalues  of 
a  complex  matrix  and  hence  it  is  of  interest  to  examine  the 
discrepancies  between  W(A)  and  the  convex  hull  of  the  spectrum 
of  A.  Towards  this  end,  the  investigator  and  two  graduate 
students  have  developed  FORTRAN  programs  for  plotting  WfA)  and 
some  of  the  higher  numerical  ranges 

n 

W(A,c)  *  f  £  c  x*Ax  I  x_,...,x  orthonormal 1 
jpij  mm  m  1  n 

where  c, , . . . ,c  are  complex  numbers.  Currently  we  are 

1  n 

examining  PASCAL  implementations  of  these  programs  which  might 
be  easily  translated  into  the  Department  of  Defense's 
forthcoming  language  ADA. 

An  interesting  result  concerning  W(A,c)  which  has  been 


.  \ 
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established  under  the  current  grant  characterizes  the  class  of 
essentially  hermitian  matrices  whenever  c,  +•••+  c  is  nonzero 
and  at  least  two  of  the  c's  are  distinct.  'This  extends  a 
result  originally  due  to  0.  Taussky  f 131 . 

Another  topic  of  research  in  this  area  concerns  the 
geometry  of  the  G-bilinear  range.  If  1.  &  r  <  n  and  if  G  ** 
[0^62:62]  is  an  r  by  3r  complex  matrix  then  W(A,G)  is  defined 
to  be  the  totality  of  sums 

y*Ax  +  •  •  •  +  y*Ax, 

where  xjf . . . ,xf ,y  j, . . . ,yr  run  over  all  n-tuples  of  complex 
numbers  satisfying 

[x’Jxj  -  Gj,  [y*yil  =  G2,  ryjx^  =  g3. 

In  addition  to  the  results  concerning  the  convexity  of  WfA,G) 
reported  in  T91 ,  the  authors  have  obtained  bounds  on  the 
maximum  modulus  of  any  element  in  W(A,G).  These  results  have 
since  appeared  in  T81. 

Results  have  also  been  obtained  concerning  the  special 
case  r=l,  G=n_,lrql,  |ql  s  1.  Several  papers  concerned  with 
this  case  have  appeared  fl,5,8r15i.  Computer  proqrams  to  plot 
approximations  to  W(AfG)  in  this  case  have  been  developed  and 
are  currently  being  expanded. 
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All  of  the  above  results  were  surveyed  at  the  SIAM 
conference  on  Applied  Linear  Alqebra  which  was  held  in  April 
of  1982. 

3.  Inequali ties  for  Principal  Suhmatri ces  of  Positive 
Def i ni te  Hermitian  Matrices  . 

In  f6J,  a  simple  proof  of  a  result  due  to  J.  Chollet  f 2 1 
is  given.  Specifically/  let  A  be  an  n-square  positive 
definite  hermitian  matrix.  Then  for  every  integer  la;  r  £  n 
and  every  integer  sequence  1  £  w^  £  •  •  •  £  w^_  £  n  the  matrix 

A_3-fwl  -  (Atwl)"1 

is  positive  semidefinite  hermitian.  Here  ATwl  denotes  the 
principal  submatrix  of  A  lying  in  the  rows  and  columns  labeled 
by  w.  This  extends  a  result  originally  due  to  I.  Schur  fill. 


4.  Generalized  Matrix  Functions  and  The  Index  of  a  Symmetry 
Class  of  Tensors  . 

This  past  period  of  research  has  also  seen  the 
publication  of  some  joint  work  with  J.  Chollet  171  in  which  we 
show  that  the  m*"*1  Exterior  product  of  an  n  m-1) 
dimensional  vector  space  is  the  only  symmetry  class  of  tensors 
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with  the  property  that  a  decomposable  tensor  x .*•••*  x  is 

1  m 

zero  when  the  dimension  of  the  linear  span  of  x, ,...,x  is 

1  m 

m-1.  This  property  of  the  exterior  product  is  then  abstracted 
to  obtain  a  definition  of  the  index  for  a  symmetry  class  of 
tensors.  The  proof  of  this  result  is  based  on  several 
interesting  characterizations  of  the  determinant  and  permanent 
functions  amoung  the  class  of  generalized  matrix  funtions 
induced  by  a  group  and  a  character. 
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ABSTRACT 

R.C.  THOMPSON:  NORMAL  MATRICES,  SINGULAR  VALUES,  SMITH 
INVARIANTS,  HANKEL  AND  TQEPLIT?.  MATRICES 


The  author's  research  during  the  period  October  1,  1981  - 
September  30,  1982  was  concerned  with:  (1)  eigenvalues  and 
diagonal  elements  of  normal  matrices,  t 2*  Smith  invariants  of 
matrices,  (3)  matrix  valued  inequalities,  and  (4)  Hankel  and 
Toeplitz  matrices.  There  was  a  one  hour  invited  address  on 
Core  Linear  Algebra  at  the  SIAM  meetinq  on  Applied  Linear 
Algebra  at  Raleigh,  April  1982,  and  a  one  half  hour  invited 
address  on  Smith  invariants  at  the  Linear  Algebra  meeting 
immediately  preceding  the  American  Math.  Society  meeting  in 
Toronto,  August  1982. 
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1.  Diagonal  elements  of  normal  matrices 

The  relationship  between  the  diagonal  elements  and  the 
eigenvalues  of  a  normal  matrix  N  is  a  long  standing  open 
question,  with  applications  in  elementary  particle  physics. 
Let  dj,...,dn  be  the  diagonal  elements  and  Aj,...,An  the 
eigenvalues.  The  following  inequality  was  proved  by  R.C. 
Thompson  some  years  ago.  For  every  complex  number  z, 

Y.  I*  -  d.  I  -2  Bin  I*  -  d.  I  £  |z  -  A  |  -  2  min  |s  -  A  |.  (l) 

i«l  Uiin  x  i.l  1  l£i£n  x 


This  is  an  infinite  class  of  conditions  relating  diagonal 
elements  and  eigenvalues.  It  is  also  clear  that 


-  s 


(2) 


where  S  is  a  doubly  stochastic  matrix  constructed  from  the 
genvectors  of  N.  It  was  conjectured  in  the  proposal,  and  in 
research  problem  that  appeared  in  Linear  and  Multilinear 
Algebra,  that  (1)  and  (2)  are  not  adequate  to  characterize  the 
diagonal  elements  of  a  normal  matrix.  This  has  now  been 
established,  in  a  paper  accepted  by  the  journal  Linear  Algebra 
and  Its  Applications. 
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2.  Smith  Invariants  of  matrix  sums 

A  paper.  Author  vs  Referee,  somewhat  expository, 
discussing  the  problem  of  how  invariant  factors  behave  when 
matrices  of  integers  add,  has  been  submitted  to  a  iournal. 

3.  Matrices  of  algebraic  integers 

Jointly  with  M.  Newman,  a  manuscript  has  been  prepared 
discussing  the  Smith  invariants  and  the  eigenvalues  of 
matrices  of  algebraic  integers.  This  sort  of  research  is 
significant  since  it  can  shed  light  on  the  behaviour  of 
singular  values.  The  manuscript  is  still  undergoing  revision. 

4.  L.  Puctyl/ni kov  published  certain  assertions,  in  two 
papers,  concerning  the  structure  of  the  spaces  of  Hankel  and 
Toeplitz  matrices,  without  proof.  These  matrices  have  many 
applications  in  signal  processing  and  in  control  theorv. 
Proofs  have  been  constructed  for  Puctyl 'nikov's  assertions 
from  his  first  paper.  Work  is  underway  concerning  the 
material  in  the  second  paper.  Along  the  same  lines,  a 
manuscript  of  150  typed  pages  has  been  prepared  translating 
and  improving  a  Russian  language  book  "Hankel  and  Toeplitz 
matrices  -  Algebraic  theory"  publ i shed  by  I.s.  Iohvidov  in 


L9 
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5.  Matrix  valued  inequalities 

Matrices  with  quaternion  elements  occur  naturally  in 


representation  theory. 


For  these  matrices  a  triangle 


inequality  has  been  found:  Given  a  matrix  A  of  quaternion 


entries,  let 


(AA*)1/2 


where 


denotes  conlugate 


transpose.  Then:  If  A,  B  have  quaternion  entries,  unitary 
matrices  U,  V  always  exist  such  that 


A  +  B  |  *  tt|a|tj*  +  v|b|  Vs 


This  will  appear  in  Linear  and  Multilinear  Algebra. 


6.  Invited  talks 

An  invited  one  hour  lecture  on  "Core  linear  algebra"  was 
given  at  the  Applied  Linear  Algebra  meeting  held  by  SIAM  (= 
Society  for  Industrial  and  Applied  Mathematics)  in  Raleigh, 
NC.,  in  April  1982.  A  one  half  hour  invited  talk  on  "Smith 
invariants  of  matrix  sums"  was  given  at  the  Linear  Algebra 
meeting  held  in  Toronto,  August  1982.  A  paper  based  on  the 
first  talk  is  being  prepared,  and  one  based  on  the  second  has 
been  submitted  to  a  journal.. 
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Papers  published 

1.  Cyclic  relations  and  the  Goldberg  coefficients  in  the 
Campbell-Baker-Hausdorf f  formula.  Proceedings  Amer.  Math. 
Soc.,  86,  1.2-14,  1982. 

2.  The  Jacobi -Gundelf inger-Frobenius-Iohvidov  rule  and  the 
Hasse  symbol.  Reports  in  Linear  Algebra,  44,  189-196, 
1982. 

3.  The  Smith  form,  the  inversion  rule  for  2x2  matrics,  and 
the  uniqueness  of  the  invariant  factors  for  finitely 
generated  modules.  Reports  in  Linear  Algebra,  44, 
1.97-201,  1982. 

4.  An  inequality  for  invariant  factors.  Proceedings  Amer. 
Math.  Soc.,  86,  9-11,  1982. 


Papers  in  press  (accepted) 

5.  Doubly  stochastic,  unitary,  unimodular,  and  complex 
orthoqonal  power  embeddings.  Acta  Scientiarium 
Matheraaticarcem  (Szeged) . 

6.  A  converse  to  the  Ostrowski-'i’ausskv  determinantal 

inequality.  Linear  Algebra  and  Applications. 

7.  A  counter  example  to  a  linear  algebra  conjecture.  Linear 
Algebra  and  Applications. 

B.  The  true  growth  rate  and  the  inflation  balancing 
principle,  Amer.  Math.  Monthly. 

9.  The  matrix  valued  triangle  inequality:  ouaternion 

version.  Linear  and  Multilinear  Algebra. 


Papers  submitted 


1.  Author  vs  Referee,  a  case  history  for  middle  class 
mathematicians,  submitted. 
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Manuscripts  still  in  preparation 

1.  p-adic  matrix  values  inequalities. 

2.  Sums  of  integral  matrices. 

3.  Hankel  and  Toeplitz  matrices  -  algebraic  theory,  150 
typed  pages. 

4.  Matrices  of  algebraic  integers  (•iointiy  with  M.  Newman). 
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ABSTRACT 

M.  NEWMAN:  EXACT  SOLUTIONS  OF  LINEAR  SYSTEMS 

The  author  prepared  a  program  using  p-adic  arithmetic  to 
find  the  exact  solution  of  an  integral  linear  system  of 
equations.  The  program  compares  favorably  in  execution  time 
with  a  standard  elimination  method  program  using  real 


arithmetic. 
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Detailed  Report: 

The  work  accomplished  this  period  was  the  preparation 
(for  a  serial  machine)  of  a  program  to  find  the  exact  solution 
of  a  linear  system  by  p-adic  arithmetic,  rather  than  residue 
arithmetic.  The  advantage  is  that  only  one  prime  modulus  is 
necessary,  and  the  numerous  local  solutions  modulo  differing 
primes  are  replaced  by  simple  matrix  by  vector  multiplications 
modulo  the  single  prime.  In  addition  the  Chinese  remainder 
theorem  is  not  required;  only  some  rather  simple  multi- 
precision  multiplications,  divisions,  and  additions  need  be 
performed.  The  time  reauired  to  find  the  exact  solution  in 
this  wav  compares  very  favorably  with  the  time  required  to 
find  an  ordinary  solution.  Some  disadvantages  are  that  the 
determinant  is  not  readily  computed  this  way,  and  that  no 
advantage  is  gained  over  ordinary  residue  arithmetic  in 
finding  the  inverse.  The  method  appears  to  have  great 
potential,  and  was  first  suggested  by  J.  Dixon. 
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Morris  Newman 
PUBLICATIONS 

October  1,  1981,  to  date. 


1.  A  radical  diophantine  eauation,  J.  Number  Theory  13 
(1981),  495-498. 

2.  Cyclotomic  units  and  Hubert's  Satz  90,  to  appear  in  Acta 
Arith. 

3.  Pseudo-similarity  and  partial  unit  regularity  (with  F.J. 
Hall,  R.E.  Hartwig,  and  i.J.  Katz),  to  appear. 

4.  Similarity  over  *  SL(n,F),  to  appear  in  Linear  and 
Multilinear  Algebra. 

5.  On  a  problem  of  H.J.  Ryser,  to  appear  in  Linear  and 
Multilinear  Algebra. 

6.  Lyapunov  revisited,  to  appear  in  Linear  Algebra  and  its 
Applications. 


/ 


Page  21 

ABSTRACT 

M.  GOLDBERG:  PROBLEMS  IN  STABILITY  ANALYSIS  OF  FINITE 
DIFFERENCE  SCHEMES  FOR  HYPERBOLIC  SYSTEMS  AND  RELATED  'TOPICS 

The  research  completed  by  M.  Goldberg  under  AFOSR  Grant 
79-0127  during  the  period  October  1981  -  September  1982, 

consists  of  the  following  topics:  (1)  Convenient  stability 
criteria  for  finite  difference  aoproxi mat ions  to  hyperbolic 
initial-boundary  value  problems.  (2)  Matrix  norms: 
submultiplicativi ty  of  C-numerical  radii  and  p-norms.  (3) 
Generalizations  of  the  Perron-Frobeni us  Theorem  and 
localization  of  the  eigenvalues  with  maximal  absolute  value. 
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The  purpose  of  this  interim  report  is  to  summarize  my  Air 
Force  sponsored  research  in  stability  analysis  of  finite 
difference  approximations  for  hyperbolic  ini tial-boundarv 
value  problems  and  related  topics,  during  the  period  October 
1,  1981  -  September  30,  1982. 


1.  Convenient  Stabilitv  Criteria  for  Difference 
Approximations  of  Hyperbolic  Initial-Boundary  Value  Problems 


In  the  past  few  years,  E.  Tadmor  and  I  fl2,13,151  have 
succeeded  in  obtaining  convenient  stabilitv  criteria  for 
difference  approximations  of  ini tial -boundary  value  problems 
asociated  with  the  hyperbolic  differential  system 


3u (x, t) /3t  =  A3u(x,t)/3x  +  Bu (x, t)  +  f  (x,t)  ,  xst  0,  t^  0. 


Here,  u(x,t)  is  the  unknown  vector,  A  a  Hermitian  matrix,  B  an 
arbitrary  matrix,  and  f(x,t)  a  given  vector. 

Our  difference  approximation  consists  of  a  well-behaved 
arbitrary  basic  scheme  (explicit  or  implicit,  dissipative  or 
unitary,  two-level  or  multi-level),  and  boundary  conditions  of 
a  rather  general  type.  Although  the  basic  scheme  is  assumed 
to  be  stable  for  the  pure  Cauchy  problem,  it  is  well  known 
that  the  introduction  of  boundary  conditions  may  destroy 
stability.  Thus,  the  auestion  is  whether  overall  stability  of 
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the  entire  approximation  for  the  initial-boundary  value 
problem  is  retained. 

Our  main  stability  criteria  in  ri2,131  for  the  above 
difference  approximations  were  given  in  terms  of  the  boundary 
conditions  and  were  essentially  independent  or'  the  basic 
scheme.  Such  scheme-independent  criteria  eliminate  the  need 
to  analyze  the  intricate  and  often  complicated  interaction 
between  the  basic  scheme  and  the  boundary  conditions;  hence 
providing  a  convenient  alternative  to  the  well  known  stability 
criterion  of  Gustafsson,  Kreiss,  and  Sundstrom  T 17 l  which  is 
usually  very  hard  to  check. 

In  the  past  two  summers,  I’admor  and  I  T151  have  expended 
our  stability  criteria  in  Q2,131  to  include  a  wider  range  of 
examples.  As  in  ri2,13],  our  stability  analysis  is  restricted 
to  the  case  where  the  outflow  boundary  conditions  are 
translatory,  i.e.,  determined  at  all  boundary  points  by  the 
same  coefficients.  Such  boundary  conditions  are  commonly  used 
in  practice  and,  in  particular,  when  the  numerical  boundary 
consists  of  a  single  point,  the  boundarv  conditions  are 
translatory  by  definition. 

Our  new  criteria  in  n*»1  depend  both  on  the  basic  scheme 
and  on  the  boundary  conditions;  yet  they  are  as  convenient  as 
our  scheme-independent  results  in  fl2,131.  In  fact,  in  many 
cases  our  new  results  imply  stability  almost  automatically  as 
stated,  for  example,  by  the  following  theorem: 


Let  the  basic  scheme  be  stable  for  the  pure  initial  value 
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problem,  and  let  the  boundary  conditions  be  explicit  and 
satisfy  the  von  Neumann  condition.  If  either  the  basic  scheme 
or  the  boundary  conditions  are  dissipative,  and  if  either  the 
basic  scheme  or  the  boundary  conditions  are  two-level,  then 
the  entire  difference  approximation  is  stahle. 

It  can  be  easily  verified  that  the  new  stability  criteria 
in  [15]  yield  all  the  examples  in  our  previous  papers  ri2,131 
and  much  more.  For  instance,  we  can  show  that  if  the  basic 
scheme  is  dissipative  and  if  the  outflow  boundary  conditions 
are  generated  by  oblique  extrapolation,  by  the  Box-scheme,  or 
by  the  right-sided  weighted  Euler  scheme,  then  overall 
stability  of  the  approximation  is  assured.  For  general  basic 
schemes  (dissipative  or  unitarv)  we  can  show  that  overall 
stability  holds  if  the  outflow  boundary  conditions  are 
determined  by  the  right-sided  explicit  or  implicit  Euler 
schemes.  We  also  showed,  for  example,  that  if  the  basic 
scheme  is  two-level  and  the  outflow  boundary  conditions  are 
generated  by  horizontal  extrapolation,  then  overall  stability 
holds  —  a  result  which  may  fail  for  multi-level  basic 
schemes.  'These  examples  and  others  incorporate  most  special 
cases  discussed  in  recent  literature  [2,3,12,13,16,17, 
19-21,27,28,30] . 

2.  Matrix  Norms:  Sub-Mult.ipl  icativitv  of  C-Numerical 


Radii  and  p-Norms 
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In  the  past  few  years,  E.G.  Straus  and  I  r5-10l  have 

continued  to  study  sub-mult iplicativi ty  properties  of  various 

norms  on  CnXn,  the  algebra  of  n  x  n  complex  matrices.  This 

topic,  which  is  strongly  related  to  most  fields  of  numerical 

analysis,  is  particularly  useful  in  stabilitv  analysis  of 

difference  schemes  for  partial  differential  equations  and 

error  analysis  for  numerical  solutions  of  linear  systems. 

In  our  work  we  considered  arbitrary  norms  on  C  „  which 

n*n 

may  or  may  not  be  (sub-1  mu! tipi icati ve.  Given  such  a  norm  K 

and  a  fixed  constant  mi  0,  then  obviously  N  =  mN  is  a  norm 

m 

too.  If  N  is  multiplicative,  i.e., 

m 

N  (AB)  *  N  <A)N  (B)  ,  A,B  £  C  v  , 
m  mm  nxn 

then  m  is  called  a  mu! tiplicativi ty  factor  for  N.  With  this 
definition  our  question  was  whether  multiplicativity  factors 
exist;  and  if  so,  how  to  find  them  in  order  to  obtain 
multiplicativity  by  very  simple  means.  A  formal  answer  to 
this  question  was  given  by  us  in  r 71 ; 

(i)  All  norms  on  CnXn  have  multiplicativitv  factors. 

(ii)  If  N  is  a  norm  on  C  then  m  5s  a  multiplicativity 

n*n 

factor  for  N  if  and  only  if 

m  i  mN  s  max{N(AB>:  N  (A)  =N  (B)  }. 

The  original  reason  for  introducing  the  idea  of 
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i 


multiplicativity  factors  was  to  investigate  the 
multiplicativity  properties  of  one  of  the  best  known  norms  on 
CnXn  •  namely  the  numerical  radius  [1,4,14,18,261 , 

r(A)  =  max{|x*Ax|:  x*x  «  1). 

We  found  [5,71  that  m  is  a  multiplicativity  factor  for  r  if 
and  only  if  m  ^  4,  independently  of  dimension.  This 

unexpected  result  is  of  interest  since  r  plays  an  important 
role  in  stability  analysis  of  difference  schemes  for 

mu! ti-space-dimensional  hyperbolic  initial  value  problems 
[14,22,23,291. 

Our  next  step  was  to  investigate  C-numerical  radii  which 
constitute  a  generalization  of  the  classical  radius  r,  defined 
by  us  in  [5]  as  follows:  For  given  matrices  A,C  €  CnXn,  the 
C-numerical  radius  of  A  is 

rc(A)  «*  max{ |  tr  (CU*AU)|  :  U  n  *  n  unitary). 

Evidently,  for  C=diag (1,0, . . . ,0) ,  rc  reduces  to  ry  thus  rc  is 
indeed  a  generalization  of  the  classical  radius. 

We  have  shown  [51  (compare  T241)  that  rc  is  a  norm  on 
Cn*n  —  and  so  has  rault-<P-1  tcati ^vity  factors  —  if  and  only  if 
C  is  not  a  scalar  matrix  and  tr  C  4  0.  Multiplicativity 
factors  for  the  above  r  were  found  in  [5-81  except  for  the 

w 

case  where  C  has  equal  eigenvalues.  The  project  was  completed 
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last  September  when  multiplicativity  factors  for  all  relevant 
C-radii  were  established  in  f91. 

Our  most  recent  effort  in  this  area  was  to  obtain 
multiplicativity  factors  for  the  well  known  -norms 
(1  <  p  *  • ) : 


{  l»i,jlP>1/P  -  *■ 


€  C 


n*n 


It  was  shown  by  Ostrowski  [251  that  these  norms  are 

multiplicative  if  and  onlv  if  1.  <  p  a£  2.  For  p  ^  2  we  have 

shown  [101  that  m  is  a  mul tipUcativitv  factor  for  |  A|p  if  and 

only  if  mi  n  1-2^p  j  thus,  <n  particular,  obtaining  the  useful 

result  that  n 1-2^p  I  A  I  is  a  multiplicative  norm  on  C  ,  . 

1  'p  r  n*n 

3.  Generali zations  of  the  Perron-Frobenius  Theorem  and 
Localization  of  Eigenvalues  with  Maximal  Absolute  Value, 


In  many  instances  one  is  interested  in  localizing  an 
eigenvalue  of  maximal  absolute  value  for  a  given  matrix.  The 
most  famous  result  in  this  vein  is  the  Perron-Frobenius 
Theorem  which  states  that  a  matrix  with  nonnegative  elements 
has  at  least  one  nonnegative  eigenvalue  of  maximal  absolute 
value. 

Last  summer,  E.G.  Straus  and  I  were  looking  for 
generalizations  of  this  celebrated  theorem  that  locate  an 
eigenvalue  of  maximal,  absolute  value  within  a  certain  angle  of 
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the  complex  plane  depending  on  the  angle  which  contains  the 
elements  of  the  matrix.  More  precisely,  let  AQ(a)  denote  the 
family  of  all  n  *  n  complex  matrices  whose  entries  are 
contained  in  a  sector 


S(<»)  =  (z:  |arg  z  j^aj  0  *  a  £  w  fixed). 


For  each  A  €  A  (a)  let  8(A)  denote  the  minimal  (nonnegative) 

n 

angle  for  which  the  sector  S(8)  contains  an  eigenvalue  of  A 
with  maximal  absolute  value.  Thus,  defining 

8n(a)  =  sup  (8(A):  A  6  An(a)} 

we  posed  the  problem  of  finding  6n(a)  as  a  function  of  a  and 
n. 

Since  the  Perron-Frobenius  Theorem  states  that  8n(0)=0,  a 
wishful  generalization  would  read  8nfa)=o.  This  unfortunately 
is  not  so,  as  shown  by  us  in  ill]  where  we  give  a  complete 
description  of  the  2x2  case  as  well  as  partial  results  for 
n  at  3. 

For  n=2,  for  example,  we  have 


P2(a) 


o  for  a  £  ir/4 
a  +  ir/2  for  w/4  <  a  £  ir /2 
ir  for  a  >  ir/2, 
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where  the  discontinuity  in  B2(a)  is  typical  for  all  n.  More 
expected  properties  of  Bn(a)  ares 
fi)  Bn(a)  *  ° 

(ii)  Bn(a)  is  «  nondecreasinq  function  of  o  and  n. 
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HENRYK  MINCs  PERMANENTS,  NONNEGATIVE  MATRICES,  AND 

HADAMARD  MATRICES 


The  author's  research  was  mainly  concerned  with 
techniques  developed  by  Egorycev  and  Faiikman  in  their 
solutions  of  the  van  der  Waerden  permanent  conjecture.  Their 
proofs  were  studied  and  the  methods  were  applied  to  extend  the 
result  to  the  problem  of  determining  the  minimum  permanent  on 
a  face  of  the  polyhedron  of  doublv  stochastic  n  x  n  matrices. 

The  author  also  studied  Frechet's  problem  of  finding  a 
simple  analytic  (or  at  least  asymptotic)  expression  for  the 
minimum  determinant  of  n  *  n  (l,-l)-matrices. 
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(I)  VAN  DER  WAERDEN  PERMANENT  CONJECTURE 

Research  in  the  theory  of  permanents  and  its  applications 
in  the  last  two  years  has  been  strongly  influenced  by  the 
recent  solutions  of  the  van  der  Waerden  permanent  conjecture 
by  Egorycev  and  by  Falikman.  They  both  proved  that 

per  (S)  i  n!/nn  (1) 

for  any  n  x  n  doubly  stochastic  matrix.  In  addition,  Egorycev 
showed  that  equality  can  hold  in  (1)  if  and  o.ily  if  S  is  the 
n  x  n  matrix  all  of  whose  entries  are  1/n.  In  paper  141  a 
version  of  Egorycev"s  proof  is  given.  Tn  [3]  two  alternative 
variants  of  his  proof  are  presented.  In  T 51  versions  of 
Egorycev" s  proof  and  of  the  more  difficult  Falikman's  proof 
are  given  in  detail. 

(II)  MINIMUM  OF  THE  PERMANENT  OF  A  DOUBLY  STOCHASTIC  MATRIX 
WITH  PRESCRIBED  ZERO  ENTRIES 

Paper  17]  is  a  study  of  properties  of  matrices  with 
minimum  permanent  in  a  face  of  &n,  the  polyhedron  of  doubly 
stochastic  n  x  n  matrices,  i.e.,  for  doubly  stochastic 
matrices  with  zero  entries  in  prescribed  fixed  positions. 
Egorycev  proved  that  all  permanental  cofactors  of  a  matrix 
with  minimal  permanent  in  n  are  equal.  This  implies  that  in 


such  matrix  any  pair  of  rows  (columns)  can  be  replaced  by 

their  mean  without  change  in  permanent.  This  averaging 

process  leads  to  a  proof  of  the  van  der  Waerden  conjecture. 

Unfortunately,  in  the  case  of  doubly  stochastic  matrices  with 

prescribed  zero  entries  such  averaging  method  has  only  a 

restricted  application  (viz.,  to  rows  (columns)  with  same 

prescribed  zero  pattern) .  In  [7]  we  show  that  permanental 

cofactors  of  a  matrix  with  minimal  permanent  in  a  face  of  8n 

cannot  exceed  the  permanent  of  the  matrix,  and  the  permanental 

cofactors  of  entries  which  are  not  prescribed  are  actually  all 

equal  to  the  permanent  of  the  matrix.  This  result  is  then 

used  to  obtain  minimum  permanents  in  faces  of  fi  in  which  all 

n 

prescribed  zeros  are  restricted  to  two  rows  or  columns,  or  in 
which  the  prescribed  zeros  form  a  submatrix. 

In  some  cases  in  which  prescribed  zeros  are  located  in 
many  rows  (columns) ,  Falikman's  method  seems  to  be  more 
appropriate  than  that  of  Egorycev.  Some  partial  results,  as 
yet  unpublished,  have  been  obtained. 

(Ill)  HADAMARD  MATRICES 

Let  a (n)  be  the  maximum  of  the  determinant  of  n  *  n 
(l,-l)-matrices.  Frechet  asked  if  there  exists  a  simple 
analytic  expression  for  A (n)  as  a  function  of  n,  and  he 
proposed  the  problem  of  determining  an  analytic  asymptotic 
expression  for  A  (n)  .  It  is  known  that  A(n)  s  n11^2  ,  and  that 
equality,  for  n  >  can  hold  only  if  n  =  0  mod  4.  In  [6] 
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relevant  known  results  for  all  n  are  surveyed,  and  it  is 
concluded  that  Frechet"s  question  should  be  answered  in  the 
negative,  although 


n/2 


log  A(n)  ~  log  n 
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